Observation of gravity-capillary wave turbulence 
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We report the observation of the cross-over between gravity and capillary wave turbulence on the 
surface of mercury. The probability density functions of the turbulent wave height are found to be 
asymmetric and thus non Gaussian. The surface wave height displays power-law spectra in both 
regimes. In the capillary region, the exponent is in fair agreement with weak turbulence theory. 
In the gravity region, it depends on the forcing parameters. This can be related to the finite size 
of the container. In addition, the scaling of those spectra with the mean energy flux is found in 
disagreement with weak turbulence theory for both regimes. 
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Wave turbulence, also known as weak turbulence, is ob- 
served in various situations: internal waves in the ocean 
surface waves on a stormy sea 0|, Alfven waves in 
astrophysical plasmas Langmuir waves Q and ion 
waves [5| in plasmas, spin waves in solids. It has been 
also emphasized that wave turbulence should play an im- 
portant role in nonlinear optics [(|. However, wave turbu- 
lence experiments are scarce. Most of them concern cap- 
illary or gravity waves. For short wavelengths, capillary 
wave turbulence has been observed by optical techniques 
0, B H E3] ■ It has been reported that the height of the 
surface displays a power-law frequency spectrum/ -17 / 6 
in agreement with weak turbulence (WT) theory ll| and 
simulations For longer wavelengths, gravity wave 

turbulence has been mainly observed in situ (i.e. on the 
sea surface or in very large tanks) with wind-generated 
waves leading to power-la w sp ectra /~ 4 Q in agreement 
with isotropic WT theory [13| and simulations [14J . How- 
ever, when the turbulence is not forced by wind or by an 
isotropic forcing, mechanisms of energy cascade in the 
inertial regime c hang e, as well as the scaling law of the 
spectrum 1J, lia, [l6| , and are still a matter of debate. 



Besides scalings with respect to frequency or wave 
number, Kolmogorov-type spectra also depend on the 
mean energy flux e cascading from injection to dissipa- 
tion. This dependence is related to the nature of non- 
linear wave interactions which are different in capillary 
(3-wave interactions) versus gravity (4- wave interactions) 
regimes [ll|, EH- To our knowledge, the mean energy flux 
has never been measured in wave turbulence and no ex- 
periment has been performed to study how spectra scale 
with e. Matching the gravity and capillary spectra and 
WT theory breakdown are other open questions 17 , 18| . 
We report in this letter how power-law spectra in the 
gravity and capillary ranges depend on the forcing pa- 
rameters of surface waves. We measure the mean energy 
flux e and show that, although the scaling of the spec- 
tra with respect to frequency looks in agreement with 



WT theory in some limits, their scaling on e differ from 
theoretical predictions. 

The experimental setup consists of a square plastic ves- 
sel, 20 cm side, filled with mercury up to a height, h 
(h = 18 mm in most experiments) (see Fig. [T]). The 
properties of the fluid are, density, p = 13.5 10 3 kg/m 3 , 
kinematic viscosity, v — 1.15 10 -7 m 2 /s and surface ten- 
sion 7 = 0.4 N/m. Contrary to the usual bulk excitation 
of waves by Faraday vibrations [7], l8| , surface waves are 
generated by the horizontal motion of two rectangular 
(10 x 3.5 cm 2 ) plunging PMMA wave makers driven by 
two electromagnetic vibration exciters (BK 4809) via a 
power supplied (Kepco Bop50-4A) . The wave makers are 
driven with random noise excitation, supplied by a func- 
tion generator (SR-DS345), and selected in a frequency 
range - fdriv with fdriv — 4 to 6 Hz by a low-pass 
filter (SR 640). This corresponds to wavelengths of sur- 
face waves larger than 4 cm. This is in contrast with 
most previous experiments on capillary wave turbulence 
driven by one excitation frequency 0, @, [T3|. Surface 
waves are generated 2.2 cm inward from two adjacent 
vessel walls and the local displacement of the fluid in re- 
sponse to these excitations is measured 7 cm away from 
the wave makers. A capacitive wire gauge, perpendicular 
to the fluid surface at rest, is made of an insulated copper 
wire, 0.1 mm in diameter. The insulation (a varnish) is 
then the dielectric of an annular capacitor with the wire 
as the inner conductor and mercury as the outer one. 
The capacitance is thus proportional to the fluid level. 
A low-cost home-made analogic multivibrator with a re- 
sponse time 0.1 ms is used as a capacitance meter in the 
range - 200 pF. The linear sensing range of the sensor 
allows wave height measurements from 10 /im up to 2 
cm with a 20 mm/V sensitivity. Although resistance or 
capacitance wire probes are widely used to get precise 
measurements of the level of quasi-static liquids, their 
dynamical response in the case of a rapidly varying wavy 
surface is not well known due to possible meniscus ef- 
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fects Thus, we have hrst checked our results with 

measurements performed with eddy current displacement 
transducers or with an optical determination of the local 
slope of the surface [20| . 

The mean energy flux injected by the wave makers 
and dissipated by viscosity is determined as follows. 
The velocity V(t) of the wave maker is measured using 
a coil placed on the top of the vibration exciter (see 
Fig- HI)- The e.m.f. generated by the moving permanent 
magnet of the vibration exciter is proportional to the 
excitation velocity. For a given excitation band width, 
the rms value, ay, of the velocity fluctuations of the 
wave maker is proportional to the driving voltage U rms 
applied to the vibration exciter. The force F(t) applied 
by the vibration exciter to the wave maker is measured 
by a piezoresistive force transducer (FGP 10 daN). 
The power injected into the fluid by the wave maker is 
I(t) = —Fn(t)V(t) where Fn(t) is the force applied by 
the fluid on the wave maker. It generally differs from 
F(t)V(t) which is measured here because of the piston 
inertia. However, their time averages are equal, thus 
(/) = (F(t)V(t)). 
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FIG. 2: (color online). Probability density functions of 
the wave-height, r), for the maximum excitation amplitude 
(Urms = 0.9 V) and for 6 different values of the fluid depth, 
from h = 18, 35, 55, 80, 110 to 140 mm (see the arrow). The 
frequency band is < / < 6 Hz. Inset: Same PDFs displayed 
using the reduced variable r\ / ■J (r/ 2 ). Gaussian fit with zero 
mean and unit standard deviation ( ). 




FIG. 1: (color online). Schematic view of the experimental 
setup showing a typical time recording of the surface wave 
height, r/(t), at a given location during 50 s. (77) ~ 0. 

A typical recording of the surface wave amplitude at 
a given location is displayed in the inset of Fig. [I] as a 
function of time. The wave amplitude is very erratic with 
a large distribution of amplitudes. The largest values of 
the amplitude are of the order of the fluid depth, whereas 
the mean value of the amplitude is close to zero. 

The probability density function (PDF) of the surface 
wave height, 77, is found to be Gaussian at low forcing 
amplitude (not shown here), whereas it becomes asym- 
metric at high enough forcing (see Fig. [5]). The positive 
rare events such as high crest waves are more probable 
than deep trough waves 21 1. This can also be directly 



observed on the temporal signal 17(f) shown in the inset 
of Fig. [TJ A similar asymetrical distribution is observed 
when using water instead of mercury, although the menis- 
cus has an opposite concavity. As shown in Fig. [2j the 
asymmetry is enlarged when the largest trough to crest 
amplitudes become comparable to the height of the layer. 
However, it persists in the limit of deep water waves. 
Note that the mean value (rf) remains close to zero and 
the PDFs of the reduced variable ^/(rj 2 ) 1 ^ 2 roughly col- 
lapse (see inset of Fig. [2]). 

The power spectrum of the surface wave amplitude is 
recorded from 4 Hz up to 200 Hz and averaged during 
2000 s. For small forcing, peaks related to the forcing 
and its harmonic are visible in the low frequency part of 
the spectrum in Fig. [3] At higher forcing, those peaks 
are smeared out and a power-law can be fitted. At higher 
frequencies, the slope of the spectrum changes, and a 
cross-over is observed near 30 Hz between two regimes. 
This corresponds to the transition from gravity to capil- 
lary wave turbulence. At still higher frequencies (greater 
than 150 Hz), viscous dissipation dominates and ends the 
energy cascade. For a narrower frequency band of exci- 
tation (0-4 Hz), similar spectra are found but with a 
broader power-law in the gravity range (see inset of Fig. 
EJ). When the two wave makers are driven with two noises 
with different band widths, e.g. - 4 Hz and 0-6 Hz, the 
harmonic peak is no longer present, and gravity spectra 
display a power-law even at low driving amplitude. 

For linear waves, the cross-over between gravity and 
capillary regimes corresponds to a wave number k of the 
order of the inverse of the capillary length l c = Wj] (fig), 
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FIG. 3: (color online). Power spectra of the surface wave 
height for two different driving voltages U rm s ~ 0.2 and 0.9 
V (from bottom to top). The frequency band is < / < 6 Hz. 
Dashed lines have slopes -4.3 and -3.2. Inset: The frequency 
band is < / < 4 Hz, and U rm a = 0.9 V. Dashed lines had 
slopes of -6.1 and -2.8. 



i.e. to a critical frequency, f c = ^/g/2l c /Tr, where g is 
the acceleration of gravity. For mercury, l c — 1.74 mm 
and f c ~ 17 Hz corresponding to a wavelength of the 
order of 1 cm. The insets of Fig. [3] and Fig. [4] show a 
correct agreement in the case of a narrow driving fre- 
quency band. We also observe that the cross-over fre- 
quency increases with the driving amplitude and with 
the width of the driving frequency band (see the inset 
of Fig. [4j. This can be due to the fact that the above 
estimate of f c is only valid for linear waves. The capil- 
lary length cannot be significantly changed using other 
interfaces between simple liquids and air. It is at an in- 
termediate scale between the size of the experiment and 
the dissipative length. In this laboratory-scale experi- 
ment, this limits both the gravity and capillary regimes 
to less than a decade in frequency. With laboratory-scale 
experiments, we can study full range gravity waves with 
a liquid- vapor interface close to its critical point and full 
range capillary waves in a micro-gravity environment. 

Surface wave turbulence is usually described as a con- 
tinuum of interacting waves governed by kinetic-like 
equations in case of small nonlinearity and weak wave 
interactions. WT theory predicts that the surface height 
spectrum S n (f), i.e. the Fourier transform of the au- 
tocorrelation function of 7/(i), is scale invariant with a 
power-law frequency dependence. Such a Kolmogorov- 
like spectrum writes 



S V (J) « ^<?r 4 



/ s for capillary waves [11], 
for gravity waves [13], 



(1) 



where e is the energy flux per unit surface and density 
[S v (f) has dimension L 2 T and e has dimension (L/T) 3 }. 
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FIG. 4: Slopes of surface-height spectra for gravity (full sym- 
bols) and capillary (open symbols) waves for different forcing 
band widths and intensities: (o) to 4 Hz, (v) to 5 Hz and 
(□) to 6 Hz. Power-law exponents of gravity wave spectrum 

( — ■) and capillary waves spectrum ( ) as predicted by WT 

theory (Eq. [1]). Inset: Cross-over frequency between gravity 
and capillary regimes as a function of the forcing intensity 
and band width. 



In both regimes, these frequency power-law exponents 
are compared in Fig. [4] with the slopes of surface height 
spectra measured for different forcing intensities and 
band widths. The experimental values of the scaling 
exponent of capillary spectra are close to the expected 
/~ 2 ' 8 scalin g as already shown with one driving fre- 
quency 0, fslTTloj] or with noise [13] ■ Figure 3] shows that 
this exponent does not depend on the amplitude and the 
frequency band of the forcing, within our experimental 
precision. For the gravity spectrum, no power-law is 
observed at small forcing since turbulence is not strong 
enough to hide the first harmonic of the forcing (see 
Fig. [3]) . At high enough forcing, the scaling exponent of 
gravity spectra is found to increase with the intensity 
and the frequency band (see Fig. SJ) . For gravity waves, 
the predicted /~ 4 scaling of Eq. (fTJ) is only observed for 
the largest forcing intensities and band width (see Fig. 
2J). The dependence of the slope of the gravity waves 
spectrum on the for cing characteristics can be ascribed 
22]. Similar results in the gravity 
found in a much larger tank 
with sinusoidal forcing [23j . 



to finite size effects 
range have been recent! 



We finally consider how those spectra scale with the 
mean energy flux e = (I)/(pSp) where (/) is the mean 
power injected by the wave maker and Sp is the area of 
the wave maker. With given ay, we have first check that 
(I) in proportional to Sp and decreases by a factor 13 
when mercury is replaced by water. Our measurements 
also show that (/) cx o v with a proportionality coeffi- 
cient of order 10 W/(m/s) 2 (see the inset of Fig. [5]). We 
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thus have e oc cay where c has the dimension of a veloc- 
ity. If we assume that e should involve only large scale 
quantities, it cannot depend on surface tension or viscos- 
ity. Then c is a characteristic gravity wave speed at large 
wave length. The dependence of e on c can be ascribed to 
finite size effects. The inverse travel time of a wave within 
the tank or the frequency difference between the discrete 
modes of the tank both scale with c. Discreteness also 
explains why a spectrum with enlarged peaks is observed 
in Fig. |4] at low forcing in the gravity range. However, 
the large enough values of e required to observe power 
laws, are more than one order of magnitude smaller than 
the critical flux (jg/p) 3 ^ 4 ~ 2200 fcm/s) 3 corresponding 
to the breakdown of weak turbulence 1181 1 . 
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FIG. 5: (color online). Spectra of the surface wave amplitude 
divided by the variance Oy of the velocity of the wave maker 
for different forcing amplitudes, ov = 2.1,2.6,3.5 and 4.1 
cm/s. The frequency band is < / < 4 Hz. The dashed line 
has slope —5.5 whereas the full line has slope —17/6. The 
mean injected power is displayed as a function of ay in the 
inset. The best fit gives a slope 11.5 W/(m/s 2 ). 

The best choice in order to collapse our experimental 
spectra on a single curve for different values of try is dis- 
played in Fig. [5] where the power spectral density divided 
by <jy is plotted versus /. Surprisingly, spectra are col- 
lapsed on both the gravity and capillary ranges by this 
single scaling. Their dependence on the mean energy flux 
e thus corresponds neither to prediction of WT theory for 
the capillary regime (e 1 / 2 ) nor to the one related to the 
gravity regime (e 1 / 3 ) but is linear with e. This discrep- 
ancy can result from several reasons. First, the size of 
the container is too small to reach a forcing-independent 
gravity regime. Second, capillary and gravity regimes 
probably interact such that it may be wrong to consider 
them independently as in Eq. ([I]). Third, we observed 
that the energy flux strongly fluctuates and takes both 
positive and negative instantaneous values much larger 
than its mean. The possible effect of these fluctuations 



on wave turbulence deserves further studies. 
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